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The field of a point electric dipole in an infinite dielectric is obtained by placing 
the dipole at the center of a spherical cavity of radius R inside the dielectric and then 
letting R — > 0. The result disagrees with the elementary answer found in textbooks. 
The mathematical and physical reasons for the disagreement are discussed. The 
discrepancy is confirmed by the same limiting procedure applied to a uniformly 
polarized sphere embedded in the dielectric. We next solve the same problem for a 
polarized spheroid immersed in an infinite dielectric and find that the asymptotic 
potential shows an unexpected shape dependence, even after taking the limit of an 
arbitrarily small spheroid. By considering both oblate and prolate spheroids and 
taking appropriate limits, we recover either the elementary textbook answer or the 
previous result found for the polarized sphere. 

I. INTRODUCTION 

Historically, electromagnetism, and particularly electrostatics, has been a rich source of 
beautiful mathematical physics problems, most of which are quite standard by now. Yet, 
from time to time, a closer look at certain simple and seemingly exhausted problems might 
surprise even the experienced practitioner. We start by discussing the elementary problem of 
determining the electrostatic field produced by a pure (point) dipole embedded in the bulk 
of an infinite linear dielectric medium. This problem is solved by two apparently equivalent 
methods. The first makes use of an elementary argument found in textbooks and the other 
consists of putting the dipole at the center of a spherical hole in the dielectric and then 
letting the radius of the hole tend to zero. The discrepancy between the results might 
surprise the reader as much as it surprised the authors. The discrepancy is corroborated 
by the same limiting procedure applied to a uniformly polarized sphere embedded in the 
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dielectric. 



Next we solve for the electrostatic field of a uniformly polarized spheroid in an infinite 
dielectric. This solution is an interesting exercise in mathematical physics involving in a 
simple way Legendre functions of the second kind, which are seldom used in the standard 
electromagnetism textbooks. We find that the asymptotic potential exhibits a shape depen- 
dence. By taking appropriate limits, we recover either the elementary textbook answer or 
the previous result found for the polarized sphere. The dependence of the electrostatic po- 
tential on the shape of the spheroid, even after taking the limit in which the spheroid shrinks 
away keeping a finite dipole moment, is unexpected and to a certain extent non- intuitive. 
This physical effect appears to have been overlooked by standard textbooks. 



The problem of obtaining the field produced by a dipole in a dielectric medium is one 
of those elementary problems that is present (solved or proposed) in a variety of textbooks. 
The well-known solution is trivial. The physical dipole consists of two opposite point charges 
(q and —q) separated by the distance d. Letting d — > with qd = p gives the pure dipole. 
Because in a linear dielectric medium Gauss' law § D • da = q establishes that each of the 
charges q and —q will be screened by polarization charges to q' = qto/e, the dipole moment 
will be screened by the same factor, so that the actual (effective) dipole moment is 



Equation (jl|) is the answer found in standard textbooks (see Ref. |l]for example). In other 
words, for a point dipole parallel to the z-axis and located at the origin, the electrostatic 
potential inside the infinite linear dielectric medium in spherical coordinates (r, 9, (p) is 



with p given by Eq. ([!]). 

Now let us solve the "same" problem by putting the pure dipole p at the center of an 
empty spherical hole of radius R cut out of the dielectric medium and then letting R — > 0. 
It is appropriate to make use of the general solution to Laplace's equation in spherical 
coordinates for problems with azimuthal symmetry. It is easy to see that the boundary 
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conditions can be satisfied by taking only the £ = 1 term of the azimuthally symmetric 
general solution, so the electrostatic potential inside the hole is 

Po cos 9 



$ (1) (r) = Ar cos 6 + 
and the potential outside is 
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Note that inside the hole the singular term corresponds to the pure dipole singularity with 
dipole moment po, because the dipole is in vacuum. Outside, only the term that decreases 
with r is present, with the factor p' to be determined. By requiring the continuity of the 
scalar potential (equivalent to the continuity of the tangential component of the electric 
field) and of the radial component of the electric displacement vector (D^ = — eoV^ 1 -* and 
= — eV^ 2 -*) at the boundary r = R, we obtain 
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The application of the boundary conditions in Eqs. (£|) and 
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leads to the equations 
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whose solution is 
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According to Eqs. (|8]) and (|j), the electrostatic potential outside the hole is that of a point 
dipole in vacuum with effective dipole moment p' given by Eq. @. In the limit R — > 0, the 
dipole potential everywhere except at the origin is given by Eq. (|3~T5| ) with p' determined by 
Eq. (Q). Surprisingly, this effective dipole moment disagrees with the one given in Eq. ([!]) 
by means of the previous elementary argument. 
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The reason for the discrepancy appears to be the lack of commutativity of two successive 
limits. The result ([I]) corresponds to putting the two opposite charges outside the hole in 
the dielectric, letting the radius of the hole tend to zero first, and then making the distance 
between the charges arbitrarily small, thus creating a dipole at the origin. To obtain the 
result (0), we first let the distance between the charges tend to zero, creating a point dipole 
at the center of the hole, and only later do we make the radius of the hole arbitrarily small. 
A physical explanation for the discrepancy is that in the first case, but not in the second 
case, the charges are always screened by the dielectric. 

We might argue that the dipole moment associated with the polarization charges on the 
surface of the hole added to po leads to a total dipole moment given by Eq. @, which is 
in fact vindicated by an explicit calculation. This argument, however, misses the point. 
The surprise comes from the fact that, if only the free dipole moment po is considered, its 
reduction by the dielectric constant factor does not account for the screening effect due to 
the polarization of the medium. This behavior contrasts sharply with that of a point charge 
at the center of the hole, whose field in the interior of the dielectric is obtained by simply 
replacing the free charge q by qeo/e in the vacuum field. 



III. UNIFORMLY POLARIZED SPHERE IN A DIELECTRIC 

To check the previous result in Eq. (0) and allow for a generalization in Sec. [IV], let us 
consider a uniformly polarized sphere (electret) of radius R, with polarization P along the 
z axis, P = -Pok, surrounded by an infinite dielectric whose dielectric constant is e. The 
potential has no singularity inside the sphere, so we have 

$W ( r ) = Br cos 6 (0 < r < R) (10a) 

for the potential inside the sphere, while the potential outside is 

$0)(r) = / CQ f . (r>R) (10b) 
47reo r A 

We now notice that = — eoV^ 1 - 1 +P and apply the same boundary conditions as before 
to obtain 

and 
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which are solved by 



and 



P - = B, (12) 



B = W (13 » 
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The resulting electrostatic potential inside the polarized sphere is 

2e + e 



$W(r) = P ° r cos 9, (0 < r < R) (15a) 



and the potential outside is 



^f2V \ -R 3 -Po COs6> / lr i\ 

$ ( r ) = 2e + e r 2 ' (r } ( } 



If we let R — > and P ~ ^ 00 i n such a way that po = (4/3)ttR 3 P remains fixed, we would 
expect to recover the point dipole po at the origin embedded in the infinite dielectric. In 
such a limit, the potential everywhere except at the origin becomes 

= o\t \ 3en 1 Po cos 6 , . 

$( 2 ) r = — — . r>0 16 
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This result coincides with the R — > limit of the previous problem of the point dipole at 
the center of an empty sphere inside the dielectric. 

Here, again, the dipole moment of the polarization charges on the spherical surface of 
the dielectric leads to the total dipole moment Thus, the field inside the dielectric is 
obtained from the vacuum field by reducing the free dipole moment by a factor that differs 
from the screening factor for a point charge. 



IV. UNIFORMLY POLARIZED SPHEROID IN A DIELECTRIC 

To put the results of Sec. |l| in a broader context, which will make possible a further 
investigation of the origin of the discrepancy encountered above, we will examine a third 
"interpolating" problem. Consider a uniformly polarized hole (electret) in the dielectric 
medium with the shape of a spheroid (an ellipsoid of revolution). 
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A. The Oblate Case 

The oblate spheroidal coordinates are defined by (see Ref. || for example) 

x = a cosh \i sin v cos <p 

y = a cosh p sin v sin ip (17) 
z = a sinh p cos v , 

with p > 0, < v < 7T, < (p < 2n, and a a positive real number. The surface of the 
spheroid is defined by /i = p , while its interior is determined by /i < /i . It is easy to see 
that the surface of the spheroid is given in cartesian coordinates by 

£ + £ + £ = i. (is) 

where X = acosh/io and Z = asinh/io, so that X > Z. The ellipsoid is oblate, that is, 
flattened along the z direction. 
In terms of the new variables 

£ = cosv (-1<£<1) 

V - ; (19) 

( = sinh /i , (0 < ( < oo) 



we can write 



with 



x = pcosip 

y = psimp (20) 

z = <C , 

p = a[(l-a(l + C 2 )] 1/2 - (21) 



The surface of the spheroid is now given by ( = Co- Laplace's equation for the potential is 
separable in these coordinates,! and its solution with rotational symmetry about the z axis, 
which is acceptable in the present physical circumstances, is 

oo 

$ (1) (£,C) =Y t Pi(0[^tPt(iC) + B l Q t (iO] (C < Co) (22a) 
inside the spheroid, and 

oo 

$ (2) (£, = E Pt(£) [QmO + DeQSO] (C > Co) (22b) 
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outside the spheroid, where Pi is the £th Legendre polynomial and Qi is the Legendre 
function of the second kind of order £. The absence of Qe(£) is necessary to guarantee the 
regularity of $ on the z axis (£ = 1). 

An inspection of Eq. ([T7|) shows that asymptotically C plays the role of a radial coordinate. 
More precisely, for large /x, we have £ ~ cos# and £ ~ r/a with r, # spherical coordinates. 
This observation strongly suggests that the terms with I = 1 alone will suffice to satisfy the 
boundary conditions, and accordingly we take 

$ (1) (£,C) = Pi(0[APi(iC) + BQi(i()] (C < Co) (23a) 



$ (2) (£,C) = Pi(0[CPi(i()+DQ 1 (i()] , (C > Co 



where 



Pi(0 = £ and g 1 «) = Ccot- 1 c-i. 

It is not difficult to show that for large C 

QiM) - 
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(23b) 
(24) 

(25) 



Therefore, the correct asymptotic behavior of $ requires that C = 0. As in the spherical 
coordinates case, it is necessary to take B = to avoid unphysical singularities. Indeed, 
the ^-component of the electric field associated with the term Pi(£)Qi(iC) is proportional 
to hJ 1 d[P 1 {^)Qx{iO]/di = a-\l - £ 2 ) 1 / 2 (£ 2 + C 2 )~ 1/2 (Ccot- 1 C - 1), which is infinite at 
£ = C = 0, that is, at the circumference p = a on the xy-plane. Thus, we try to satisfy the 
boundary conditions with (the imaginary unit has been absorbed into the coefficient A) 



(C < Co) 



and 



$( 2 )(£,C) = J D£(Ccor 1 C-i). (C > Co) 

The continuity of the potential at the surface of the spheroid yields 

ACo = -D(Cocot- 1 Co-l). 



(26a) 



(26b) 



(27) 



The continuity of the normal component of D on the surface of the spheroid demands that 
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with 

-e+c 2 v< 2 

i + C 

Because the unit outward normal vector to the surface of the spheroid is 



\ = «(T7^ " ( 29 ) 



e 



C 



dr/d( 1 dr 



\dr/d(\ h^dC 
it follows that 



The substitution of Eqs. fl26|) and ( [31] ) into Eq. ( p8|) leads to 



-e A + P a = ~eD (^cot- 1 Co - I • ( 32 ) 

The solution for A and -D is 



(30) 



P • e ? = P k • e ? = P lk • % = pS . (31) 



A = (cot^Co-l/CoW f33l 
(eo-eJoot-iCo + ^-g' 

P = p- ? . (34) 

(e - e ) cot-% + rfe - g 

The volume of the spheroid is 

4 4 
V = -7r(acosh/i ) 2 asinh/z = -7ra 3 (l + C 2 )Co , (35) 

so that the dipole moment of the spheroid is 

Po = ^ra 3 (l + Co 2 )CoiV (36) 



For large £, we have £ ~ cos# and ( ~ r/a, so that, with the use of Eqs. (p5|) and (p6|), 
we find that the asymptotic behavior of the potential is 

$(obw n _> fo p^cos# 

(6-eo)Co(l + Co 2 )cot- 1 Co + eo + (eo-e)Co 2 47r6 r» ' 1 J 

This result is quite unexpected! Our intuition leads us to believe that seen from far away, 
it is impossible to tell a uniformly polarized sphere from a uniformly polarized ellipsoid. The 
shape independence of the asymptotic potential prevails only for the vacuum (e = e ). For 
a dielectric, the asymptotic potential depends on the shape of the dipole distribution near 
the origin. 
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B. The Prolate Case 

The analysis of the prolate case runs along similar lines. The prolate spheroidal coordi- 
nates are defined byi 

x = a sinh \x sin v cos (p 

y = a sinh ji sin v sin ip (38) 
z = a cosh n cos v . 

The surface of the spheroid is defined by p = /i , while its interior is determined by p, < p Q . 
The surface of the spheroid is again given in cartesian coordinates by Eq. ( |T8"D with X = 
a sinh fi Q and Z = a cosh fi Q , so that Z > X and the spheroid is elongated in the z direction. 
In terms of the new variables 

£ = cosw (-1<£<1) 

, \ (39) 

7] = cosh jl , (1 < T] < oo) 



we can write 



with 



x = p cos ip 

y = p sin Lf (40) 
z = a£r), 

p = 4(i-£W-i)] 1/2 - (4i) 



The surface of the spheroid is now given by r\ = r] , and the solution to Laplace's equation 
that suits our problem is 

$M{Z,r t )=FP 1 {Z)P 1 (7i), (v< Vo ) (42) 

and 

& 2 \t, V )=GP 1 (Z)Qi(v)- iv>Vo). (43) 

Here 

Qi(r/) = |ln^-l, (44) 



whose asymptotic behavior for large r\ is 



^-^^ (45) 
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The application of the boundary conditions at the surface of the uniformly polarized 
prolate spheroid yields 



1 2 r? -l VO J u 

ep-e J n VO+l _i_ eyo _ eo. 

2 J7o— 1 r$-l 170 



(46) 



G= . (47) 

„ _1 " 2 - 



2 170 — I ??o -1 W 

The volume of the prolate spheroid is 



4 4 
V" = -7r(asinh/xo) 2 acosh/i = -^(^(Vo ~~ ^)Vo , (48) 



and its dipole moment is 

Po = \™ i {r,l- 1)77 Po- (49) 

As for the oblate case, we have for large rj, £ ~ cos# and 77 ~ r/a, so that the asymptotic 
behavior of the potential is 

$(Pr) — ► ^ _Po_ costf ( } 

Once again the asymptotic potential exhibits a surprising shape dependence that is absent 
only in the case of the vacuum. For a dielectric the asymptotic potential allows us to tell the 
difference between a uniform dipole density distributed within a sphere, an oblate spheroid, 
or a prolate spheroid. If a hole in the dielectric is filled with a uniform charge density, no 
such shape dependence is observed. The total induced charge depends only on the internal 
free charge and on the dielectric constant. 

Note also the striking result that the electric field is uniform inside the spheroid because 
the potential is of the form $W = Az in both cases, as Eqs. ( |26aj) and ( ^2| ) show. 

V. LIMITING CASES AND CONCLUSION 

Let us define the screening factor a as the coefficient that multiplies the vacuum asymp- 
totic dipole field to give the asymptotic dipole field in the presence of the dielectric medium. 
From Eq. (|50|) it follows immediately that 

-1 



a = e 



—5— VoiVo ~ 1) ln 7 + e + (e - e )r7 



(51) 



11 



The spherical limit is reached by letting rj = cosh fi — > oo and a — > in such a way 
that a cosh fi = R remains fixed. Then we obtain X = Z = R in Eq. ([T^) and the spheroid 
degenerates into a sphere. If we take into account that for large ryo 

Vo + l = A + JL 

770 - 1 r?o 3»/j} 

it is easy to show that 



i»r Lr r = ^ + /s + ---> (52) 



lim a = — . (53) 

Thus our previous result for the uniformly polarized sphere is recovered. 

Let us now examine the line dipole limit, reached by letting 7] — > 1 or, equivalently, 
Hq — > 0. In this limit we have X = and Z = a, so that the ellipsoid ([18]) reduces to a line 
segment (a rod) along the z-axis. If we recall that lim^o % m % — 0, we can readily show 
that 

lim a — — , (54) 
rio-*l e 

and the standard answer (HD is regained. 

For the oblate spheroid we can also consider two limiting cases with the help of Eq. (|37|). If 
we let Co — » 00 and a — > with a^o = R, the spherical limit is reached and a — > 3e / (2e + e ), 
as it should. If Co — ^ 0, the spheroid becomes a "pancake" describing a dipole layer, and 
a — > 1. This result appears to be of some interest, inasmuch as the presence of the dielectric 
does not change the vacuum field. 

We believe that further discussion of the physical grounds for the discrepancies is neces- 
sary. The screening factor a is a measure of the total dipole moment in the presence of the 
dieletric relative to the vacuum dipole moment pq. Thus, the dipole moment of the charges 
induced on the surface of the hole depends not only on the free dipole moment po and the 
dielectric constant, but also on the shape of the hole. That the shape dependence persists 
in the infinitely small hole limit seems to be related to the singularity of the dipole field, 
which is stronger than that of the monopole field. For a shrinking finite charge distribution 
(monopole), Gauss' law forbids this effect. The reader might want to generalize the dipole 
result for higher multipole moments. 

It is suspected that such a shape dependence would manifest itself in the dynamical case, 
that is, in the radiation from a point dipole embedded in an infinite dielectric. Such a 
phenomenon might be of relevance in condensed matter physics. For instance, it might give 
rise to classical effects in the theory of quantum dots. 
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